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Abstract 

There are problems with defining the thermodynamic hmit of systems with long-range interac- 
tions; as a result, the thermodynamic behavior of these types of systems is anomalous. In the 
present work, we review some concepts from both extensive and nonextensive thermodynamic per- 
spectives. We use a model, whose Hamiltonian takes into account spins ferromagnetically coupled 
in a chain via a power law that decays at large interparticle distance r as 1/r" for a > 0. Here, we 
review old nonextensive scaling. In addition, we propose a new Hamiltonian scaled by ~^ 
that explicitly includes symmetry of the lattice and dependence on the size, N, of the system. The 
new approach enabled us to improve upon previous results. A numerical test is conducted through 
Monte Carlo simulations. In the model, periodic boundary conditions are adopted to eliminate 
surface effects. 



I. INTRODUCTION 



A good description of magnetic ordering phases illuminates concepts about the critical 
behavior and possible applications of magnetic devices. We know that the magnetic behavior 
of systems decreases when the dimensionality of a physical system decreases. This led 
to incorrect beliefs that prompted specialists to lose interest in one-dimensional systems. 
However, several theoretical and practical aspects, which appeared in the physical properties 
of one- dimensional systems, caused specialists to reverse those beliefs. 
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Recently, ferromagnetism in one dimension has been reported in various systems. Transi- 
tions between two different magnetic ordering phases are obtained using different approaches 
and considerations, and we will discuss some of these. First, microscopic anomalies lead to 
important modifications in the thermodynamic properties of systems. As a consequence, 
it is suggested that anisotropy barriers contribute to the novel effect. Second, the Gibbs 
potential of a one- dimensional metal at constant magnetization is calculated to second order 
in the screened electron-electron interaction. At zero temperature, a possible paramagnetic- 
ferromagnetic quantum phase transition was found in one-dimensional metals, which must 
be first order, js] Finally, a special Hamiltonian considers long-range interactions through a 
power law that decays at large interparticle distances. It has been shown that if the range 
of interactions decreases, the critical temperature trend disappears, but if the range of in- 
teractions increases, the trend of the critical temperature approaches the mean field approx- 
imation. The crossover between these two limiting situations is preliminarily discussed. jlo| . 
That crossover is a consequence of long-range interactions, and we use it to illustrate the 
nonextensive perspective of the thermodynamics. 

To describe the behavior of systems with long range interactions, some scalings ap- 



proaches 



12| have been introduced in the literature. In the present work, we carefully 



discuss the thermodynamic behavior of systems with microscopic long-range interactions, 
taking into account a method of implementing periodic boundary conditions (section IIIip 
via a scaling Hamiltonian (section lIVp . Previously, the critical temperature between two 
states of different magnetic ordering was obtained in a spin-^ Ising linear chain, where the 
range of interactions is, at least, comparable to the size of the chain. However, in order to 
obtain suitable thermodynamic behavior, in accordance with impositions from extensivity. 
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we improve the so called nonextensive scaling (section lIVp for a symmetric one-dimensional 
lattice. 

This paper is organized as follows: In section [III we introduce the nonextensive perspec- 
tive of thermodynamics. In section IIIIl we explain the method of implementing periodic 
boundary conditions to eliminate surface effects. In section \IV\ we review the Tsallis scaling 
to get a formalism of the thermodynamics from a scaling Hamiltonian and then explain some 
important results about the presented formalism. In section |V] we summarize our work and 
make some concluding remarks. 



II. NONEXTENSIVE THERMODYNAMICS 



In this section, we introduce some fundamental facts about the nonextensive perspective 

of thermodynamics, which can be easily illustrated by long-range interaction systems. 
A. Tsallis Scaling: 

This method is useful for scaling thermodynamic quantities that depend strongly on the 
size of the system with long-range microscopic interactions. The explicit form of Tsallis scal- 
ing appears by evaluating the internal energy per particle. We take into account interparticle 
potentials f (r) with an attractive tail that decays as: 



v(r) = — . 



(1) 



The original Tsallis scaling, from the asymptotic trend of U^, for systems in one dimension, 
is given by: 



(Ari-°)/(l - a) for < a< 1 
N* = I \og{N) for a = l (2) 

l/{a - 1) for a > 1. 

The present scaling has been revised by several authors and applied to different physical 
situations. The scaling was expressed in a most appropriate form for discrete systems llf| 
and was previously written as 2°'N*. Indeed, in the present work we reviewed this kind of 
scahng and included it best approach. 
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B. Mean Field Approximation 



We consider a one dimensional model of N spins- 1. If the coupling is long-range, surface 
effects appear and become important. This fact requires working with very large systems. 
However, surface effects can be ignored, even in moderately sized systems, by applying 
periodic boundary conditions. This is a nontrivial task when long range interactions are 
considered. However, it is possible to do in a suitable form, which is introduced in section [XTTl 
The model is described by the Hamiltonian: 

n 

H = -J2 J{i- 3)siSj, (3) 

where n is the number of particles in every cell, and = ±1 Vz. The difference {i — j) 
represents the distance between two sites. A non external magnetic field is considered. The 
coupling decays as a power law given by: 

^(^-^■) = ^, (4) 
where J is a positive parameter which measures the strength of the coupling. In previous 
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works, these kinds of systems have been discussed by several authors, jl, 2.l3l.l4l.l5l.l6l.l7l.llO, 
lisl A common behavior has been conjectured [3| for generic systems with arbitrary (long or 
short)-range interactions. With the aim of calculating the mean value of the Hamiltonian, 
we proceed in the same way as in the Eq.^. (H) = E^Ii E{N)si, where E{N) = |JA^*(s); 
the factor | ensures that we do not count the same pair of spin twice and in this view; A^* 
(defined in the Eg. (1121) ) represents the effective number of the nearest neighbor spins. The 
quantity (s) is the average spin per site. This assumption helps us recover all the results 
for systems with long-range microscopic interactions in a similar manner to the traditional 
mean field approximation. This treatment can be extended even if the external magnetic 
field is nonzero. 

So, the average spin per site of the lattice is given by (s) = tanh ^PN*J{s) . In the 
absence of an external magnetic field, the magnetization is zero for a high temperature para- 
magnetic phase, and it will be nonzero at lower temperatures where the spins have sponta- 
neously aligned. For the present system, the internal energy is given by U{N) = ^NN*J{s)'^. 
Hence, the internal energy is given by U{N) = NN*Ui. The extensive property imposes 
observables to be a linear homogeneous function of N and U (XN) = XU{N). However, when 
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long range interactions are included, this property is violated, and it is easy to show that 
thermodynamic functions are homogeneous of degree 1 + |a; — 1| for a < 1 (and logarithmic 
of AT for a = 1). 

In addition, we expect that the internal energy of a magnetic system with long-range 
interactions adopts the following form: U{S,M,N) = NN*Ui{S/N,M/N), where Ui is a 
function per particle of the entropy, S, and magnetization, M. Generally, we could write 
other extensive variables like X = V, A, L, P, ■ ■ ■ (representing the volume, area, length, and 
polarization, respectively) as follows: 

U{XS, XX, XN) = A^+l"-^lt/(S, X, N). (5) 

Hence, a problem with defining the thermodynamic limit persists for a < 1, in accordance 
with Eq.([5]). The linearity is only recovered for a > 1. As previously observed, a strong 
dependence on the size of the system obstructs the behavior of the thermodynamic functions 
and relations. 



C. Tsallis Conjecture: 

As a possible way to solve this problem, Tsallis conjectured that quantities like energy 
(like internal energy, Helmholtz and Gibbs energies, and thermodynamic potentials per 
particles) and intensive variables (like T as in temperature, H as in magnetic field) scale 
with A^*. Consistency of this conjecture is shown as follows 

Gn _ Un _ ^^N_ _ (f.\ 

N*N ~ N*N ~ iV* aT ~ iV* iV 

in the thermodynamic limit. Observables per particle and intensive variables are divergent 
when interactions are long-range. However, scaling quantities are convergent anywhere. 
This kind of scaling presents a standard thermodynamic structure because it preserves the 
Euler and Gibbs-Duhem relations [l^. By the same reasoning, it is possible to define scaling 
quantities as follows: = Gn/N*, U*^ = Un/N\ T* = T/N* and H* = H/N*. With 
these definitions the previous equation is given by 

9jL = UjL -T*— - H*— f7l 
N N N N' ^ ^ 

At this stage, we focus on some problems related to the interpretation of the Eq.(I7j). For 
instance, the measured temperature is not T*, it is T; the measured internal energy per 
particle is not [/^, it is U^, etc. 
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III. PERIODIC BOUNDARY CONDITIONS 



If the range of microscopic interactions is smaller than the size of the system, thermo- 
dynamic properties are obtained from standard calculations. However, if interactions are 
long-range, surface effects appear and begin to be important for all finite sizes of systems. 
The calculation of thermodynamic quantities must be carefully done. Long-range interac- 
tions are often defined as those that do not fall faster than where D is the space 
dimension of the system. A first approach requires increasing tremendously the size of the 
box. In general, periodic boundary conditions are applied in order to eliminate surface effects 
in the calculation of the thermodynamics properties of systems. This can be exemplified by 
a central cell, which is repeated throughout space to construct an infinite lattice. During 
the course of the present study, if a particle moves in the central cell, its periodic images 
move with the same orientation in everyone of the other cells. Thus, as a particle leaves 
the central cell, one of its images will enter through the opposite face. There are no walls 
at the boundary of the central cell, and the system has no surface. In general, particles 
interact with a potential of the following form v{z) = lim^^^oo Sfc=-L fi'd^ + ^1); where the 
summation over k represents all contributions over replicated images. Thus, summations 
can be written in the following manner: 

i: 9{\k + z\) = j:g{k-z)+g{\z\) + j:g{k + z). (8) 

k=-L k=l k=l 

Two particular cases (the 1/z and the logarithmic potential) have been discussed in a pre- 
vious workfisl. 

Certainly, forces are obtained from f{z) = —dv{z)/dz. Due to the symmetry of the 
lattice, it is easy to show that contributions to the net force on every particle from all their 
images vanish. Thus, periodic summations of forces do not practically depend on the nature 
of interactions. In this manner, the net force (/tv = —dg{z)/dz — Fl) converges as quickly 
as the sum Fl converges 



L 

E 



dg{k + z) dg{k — z 



+ 



k=i I dz ' dz \' ^ ^ 

where — 1/2 < 2; < 1/2 and k ^ 1. dg{z)/dz represents the force of a couple of particles 
in the central cell, and Fl contributions of all replicated cells in all space. If we take into 
account a potential like Eq. ([1]) , net force converges in the same manner as summations from 
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FIG. 1: Linearity of magnetization is shown for L = 0, 1 and 10 and several values of n with 
N = {2L + l)n. We plot M{Xn) versus M{n), where A = (2L + 1). We see that the slope is A. In 
the inset we plot the same in log-log scale. The slope is 1. This fact emphasizes that magnetization 
is a linear homogeneous function of the size of the system. 



Eq.®: 



-E 

k=l 



(k + z) 



(k-z) 



a+l 



k=l 



(10) 



oc 



4z (l - L-^- 



If L > 1, contributions of particles and their images are included. 
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IV. SCALING HAMILTONIAN AND NUMERICAL TEST 



We consider a symmetric chain of n spin-| in a central cell that it is replicated L times 
over both sides. So, the total number of particles = (2L + l)n. Thus, we study the 
system via the Hamiltonian of the Eq.Qand apply periodic boundary conditions in the 
manner introduced in section IIIII We write the coupling as follows 

J(K) = -y , ^ ■ (11) 

Hence, in order to obtain the scaling, we can determine it using the tendency of internal 
energy, similar to Tsallis [3]. Taking into account the symmetry of the chain and the 
continuous limit, we derive the following integral: 

Un oc 2 drg{r)v{r) 



oc 2W^)"° -\ (12) 

1 — a 



where g{r) (the pair distribution function) approaches 1 for r ^ 1. Indeed, from Eq.ffT? 
the scaling coupling can be written in the approach: 



J/2"A^^-" if a < 1 

J/2 In ifa = l (13) 

(1 -«)J/2 if a > 1 



where J^i^) measures the strength of the coupling that depends on the size n of the system. 
If we combine Eq. ffTTj) and f|T3l) and substitute the result into Eq.(l3]), then we obtain the 
scaling Hamiltonian that we used to carry out numerical results with the Monte Carlo 
procedure. At this stage, we emphasize that the scaling, in terms of A^*, is 2'^N* for a < 1 
and 2A^* for a > 1. It has been well expressed in Eq.f llSp for symmetric lattices. Previously, 



(see for instance Ref . [l l| ) the scaling coupling was defined in the old form (1 — a)J/2'^ and 
not as the new form (1 — a) J/2 (from Eq. f|T3l) ) for a > 1, for a discrete lattice. Hence, if the 
old scaling is applied to systems whose most important geometrical property is symmetry, 
then it would poorly represent the trends in thermodynamic quantities. 

Thermodynamics describes the behavior of systems with many degrees of freedom af- 
ter they have reached a state of thermal equilibrium. Furthermore, their thermodynamic 
state can be specified in terms of a few parameters called state variables. At equilibrium, 
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this method of scahng the Hamiltonian allows observables to be linear with the number of 
particles. 

A. Linearity of extensive quantities 

It is well known that thermodynamic quantities have to behave linearly relative to N, 
the size of the system. For the present specific model of a chain of spins with ferromagnetic 
long-range interactions, due to the form of the scaling Hamiltonian, we expect that the 
linearity of internal energy was completely satisfied and U{\n) = XU{n). This goes against 
the nonextensive view of the thermodynamics, which predicts a nonlinear homogeneous 
behavior of internal energy, like Eq.Q. In addition, according to the Tsallis conjecture, we 
expect that other extensive quantities also become linear with A^. In Fig. [H the linearity 
of magnetization was tested from simulations for L = 0, 1 and 10 and several values of n, 
where N = (2L + l)n. We depict several values of magnetization M{Xn) versus M{n), where 
A = (2L + 1). A simple inspection reveals linearity, i.e.M{Xn) = XM{n). Simulation points 
depict right lines that differ only in their slopes. This fact is highlighted in the inset of Fig. [1] 
in log-log scale, where we observe that all lines are parallel, whose slope is 1, an indicator of 
linearity. So, we see that the linearity of an observable is satisfied when it is different from 
the internal energy, similar to the case of magnetization. 

B. Improvement of the critical temperature 

Phase transition characterization is obtained in several ways. One of the most adequate 
approaches is to define the critical point for finite systems by the Binder method. To do 
so, we need to assess the Binder cumulants of fourth order, which are defined as u„ = 
1 - (s^)/3(s2)2. Cumulants function of the temperature, intersect at a common 

point for several sizes, n, of the system. This point is the critical temperature that depends 
on the range of interactions a. We test the trend of the critical temperature and plot it as a 
function of a using the scaling Hamiltonian, and we compare the present result to the trend 
obtained previously from the old scaling. In the present study, we carry out on a linear 
chain, where the number of particles in the central cell is 10^ < n < 10^. On the one hand, 
the effective number of particles of the samples = (2L + l)n depends on the number, L, of 
replicas. On the other hand, the computation time depends just on n, not on L. We see that 
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FIG. 2: The trend of the critical temperature as a function of a, the range of interaction. For 
< a < 1, the critical temperature coincides with the mean field value. The new trend is compared 
to the old one. In the inset we show the critical temperature as a function of 1/a 

L = 10^ is good enough, to compute in the thermodynamic limit and, eliminate unnecessary 
sources of numerical fluctuations. For < a < 1, the critical temperature approaches the 
mean field value {e.g., for a = 0.5, T^/T^'^^ = 0.985 ± 0.005). But when a > 1, results from 
different scalings do not coincide and they are different from the values predicted by the 
mean field approximation. In Fig. [2l the trend of the critical temperature is depicted as a 
function of a. We include results in the range < a < 2 for the critical temperature scaled 
in the manner presented. 
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V. SUMMARY AND CONCLUDING REMARKS 



In the mean field approximation, the state of magnetic ordering of a chain of spins 
with microscopic long-range interactions constitutes an example of ferromagnetism in one 
dimension. Because interactions are short range {e.g., first neighbors) in the standard Ising 
model, no magnetic ordering is observed in one dimension. These approaches define two 
limiting cases. In the present work, the main goal is to describe the crossover between these 
two limiting cases in a suitable manner. First, if we take a scaling Hamiltonian, the nice 
extensive property is recovered, and the thermodynamic quantities are linear homogeneous 
functions against the homogeneous function of degree 1 + | o; + 1 1 , which are obtained if the 
Hamiltonian is not scaled. Second, an improvement to the nonextensive scaling is proposed 
here to obtain a more suitable method of describing the thermodynamic behavior trend in 
this kind of systems for a > 1. In the N ^ oo limit, we compare old scaling: 2"A^*, versus 
new scaling: 2N*. Third, periodic boundary conditions can be used in the approach, using 
infinite replications of a central cell and considering the contribution over all space. It has 
been shown that forces converge rapidly. However, potentials increase with the size of the 
system. The thermodynamic limit was reached in a numerical approach, with a few particles 
in the central cell and a finite number, L, of replications. Finally, we focus our attention 
on the Tsallis scaling (not the Tsallis conjecture); which is sufficient to solve the problem of 
the loss of linearity for thermodynamic quantities like energy and intensive variables, if we 
use it conveniently in the Hamiltonian. 
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